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Abstract.  A  thermal  parametric  model  has  been  developed  for  analyzing  observed 
regional  sea  temperature  profiles  based  on  a  layered  structure  of  temperature  fields 
(mixed  layer,  thermocline,  and  deep  layers).  It  contains  three  major  components:  (1)  a 
first-guess  parametric  model.  (2)  high-resolution  profiles  interpolated  from  observed 
profiles,  and  (3)  fitting  of  high-resolution  profiles  to  the  parametric  model.  The  output  of 
this  parametric  model  is  a  set  of  major  characteristics  of  each  profile:  sea  surface 
temperature,  mixed-layer  depth,  thermocline  depth,  thermocline  temperature  gradient, 
and  deep  layer  stratification.  Analyzing  nearly  15.000  Yellow  Sea  historical  (1950-1988) 
temperature  profiles  (conductivity-temperature-depth  station,  4825:  expendable 
bathythermograph.  3213:  bathythermograph.  6965)  from  the  Naval  Oceanographic  Offices 
Master  Oceanographic  Observation  Data  Set  by  this  parametric  model,  the  Yellow  Sea 
thermal  field  reveals  dual  structure:  one  layer  (vertically  uniform)  during  winter  and 
multilayer  (mixed  layer,  thermocline.  sublayer)  during  summer.  Strong  seasonal  variations 
were  also  found  in  mixed-layer  depth,  thermocline  depth,  and  thermocline  strength. 


1.  Introduction 

The  Yellow  Sea  is  a  semi-enclosed  basin  covering  roughly 
295,000  km2  and  is  one  of  the  most  developed  continental  shelf 
areas  in  the  world  seas  [Yarutgi  and  Takahashi.  19931-  While 
the  Yellow  Sea  covers  a  relatively  large  area,  it  is  quite  shallow, 
reaching  a  maximum  depth  of  about  140  m  (Figure  I).  The 
water  depth  over  most  of  the  area  is  less  than  50  m.  The 
deepest  water  is  confined  to  a  north-south  oriented  treneh 
which  runs  from  the  northern  boundary  south  to  the  100-m 
isobath,  where  it  fans  out  onto  the  continental  break.  The 
gradients  in  slope  across  the  bottom  arc  very  small.  Such  a 
broad  and  shallow  continental  shelf  leads  to  the  fact  that  the 
water  is  readily  affected  by  seasonally  varying  atmospheric  con¬ 
ditions  such  as  heating,  cooling,  and  wind  stress.  Therefore  the 
seasonal  variation  of  the  water  masses  is  remarkably  large. 
Another  feature  of  the  depth  distribution  is  the  cast'west  asym¬ 
metry.  Extensive  shoals  (<20  m)  are  located  in  the  western 
Yellow  Sea  along  the  Chinese  coast  and  are  not  generally 
found  in  the  South  Korea  coastal  regions.  .Also,  the  50-m 
isobath  is  located  more  than  100  km  from  the  Chinese  coast 
but  onlv  about  50  km  from  the  South  Korea  coast.  This  asym¬ 
metry  in  bottom  depth  is  important  for  the  shoaling  mixed- 
layer  depth.  Furthermore,  the  hydrographic  character  of  water 
masses  in  the  Yellow  Sea  also  depends  on  the  degree  of  mixing 
of  fresh  water  originating  from  the  China  continent  river  runoff 
with  the  intrusion  of  East  China  Sea  and  Kuroshio  waters. 

Before  we  sketch  vertical  structures  of  the  Yellow  Sea  shelf 
water  masses  from  the  profiles,  we  should  ask  ourselves  how  can 
the  thermal  fields  be  described  in  terms  of  a  set  of  characteristic 
parameters?  In  this  study  we  developed  a  thermal  parametric 
model  to  analyze  15.000  historical  temperature  profiles  from 
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the  Naval  Oceanographic  Office  (NAVOCEANO)'s  Master 
Oceanographic  Observation  Data  Set  (MOODS)  from  the 
years  1950  to  1988  and  to  transform  each  temperature  profile 
into  a  set  of  characteristic  parameters  such  as  sea  surface 
temperature  (SST).  mixed-layer  depth  (MLD).  thermocline 
depth,  temperature  difference  across  the  thermocline,  and  the 
deep  layer  stratification. 

2.  Master  Oceanographic  Observation  Data  Set 
(MOODS) 

The  MOODS  is  a  compilation  of  observed  ocean  data  world¬ 
wide  consisting  of  (1)  temperature-only  profiles.  (2)  both  tem¬ 
perature  and  salinity  profiles.  (3)  sound-speed  profiles,  and  (4) 
surface  temperatures  from  drifting  buoys.  These  measure¬ 
ments  are.  in  general,  irregular  in  time  and  space.  In  this  study 
we  analyze  temperature  profiles  measured  from  a  variety  of 
instruments.  Due  to  the  shear  size  (more  than  six  million 
profiles)  and  enormous  influx  of  data  to  NAVOCEANO  from 
various  sources,  quality  control  is  a  difficult  task.  Our  study 
domain  includes  the  area  325-41“N  and  118c-i27°E;  the  data 
set  within  this  region  consisted  of  nearly  15.000  profiles  after 
rejecting  certain  data  during  quality  control.  These  primary 
editing  procedures  included  removal  of  profiles  with  obviously 
erroneous  location,  profiles  with  large  spikes,  and  profiles  dis¬ 
playing  features  that  do  not  match  the  characteristics  of  sur¬ 
rounding  profiles.  In  shallow  water  this  procedure  can  be  par¬ 
tially  automated  but  also  involves  subjective  interpretation 
because  of  the  undersampling  of  MOODS  compared  to  the 
spatial  and  temporal  variability  of  the  oceanography. 

3.  Seasonal  Variation  of  the  Atmospheric 
Forcing 

The  Asian  monsoon  strongly  affects  the  Yellow  Sea  thermal 
structure.  During  the  winter  monsoon  season  a  very  cold 
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Figure  1.  The  Yellow  Sea  bathymetry.  The  data  were  ob¬ 
tained  from  the  Naval  Oceanographic  Office  DBDB5  world 
bathymetry  database.  Numbers  show  the  depth  (meters). 

northwest  wind  blows  over  the  Yellow  Sea  as  a  result  of  the 
Siberian  high-pressure  system.  The  jet  stream  is  positioned  to 
the  south  of  the  Yellow  Sea  and  the  polar  front  to  the  north  of 
Philippines.  The  mean  surface  wind  speed  over  the  Yellow  Sea 
in  January  is  nearly  f>  m/s  (Figure  2a).  The  sea  surface  tem¬ 
perature  is  6°C  at  the  northern  extent  and  I ()°C  at  the  south¬ 
eastern  extent.  The  January  surface  air  temperature  (SAT) 
varies  from  (E  to  8°C  in  the  Yellow  Sea.  roughly  2°-6°C  cooler 
than  SST.  The  Yellow  Sea  surface  loses  heat  to  the  atmo¬ 


sphere.  The  upward  buoyancy  flux  tit  the  air-ocean  interfaee 
(thermal  forcing),  together  with  the  strong  wind  stress  (me¬ 
chanical  forcing),  generates  turbulence  and  mixes  the  surface 
water  with  the  deeper  water.  The  mixed  layer  is  at  its  deepest 
(usually  tills  the  whole  water  column)  during  winter  owing  to 
both  convection  and  wind  mixing  by  the  strong  northeast  mon¬ 
soon  winds.  The  change  of  season  begins  in  March  when  the 
surface  air  temperature  is  5°C  warmer  than  in  February.  Rapid 
weakening  of  the  Siberian  high  progresses  into  April. 

In  late  April  the  polar  front  has  moved  northward  toward 
Korea  with  warm,  moist  air  following  behind.  Numerous  fron¬ 
tally  generated  events  occur,  making  late  April  and  May  highly 
variable  in  terms  of  w  inds  and  cloud  amount.  By  May  the  daily 
high  surface  air  temperatures  rise  to  150-16CC.  During  this 
period,  storms  originating  in  Mongolia  may  cause  strong,  warm 
westerlies  carrying  yellow  desert  sand  (termed  the  "Yellow- 
Wind").  By  late  May  and  early  June  the  summer  surface  at¬ 
mospheric  low-pressure  system  begins  to  form  over  Asia.  Ini¬ 
tially.  this  low-pressure  system  is  centered  north  of  the  Yellow- 
Sea.  producing  westerly  winds.  In  late  June  this  low  begins  to 
migrate  to  the  west,  setting  up  the  southwest  monsoon  that 
dominates  the  summer  months.  The  winds  remain  variable 
through  June  until  strengthening  of  the  Manchurian  low- 
pressure  system  occurs.  The  jet  stream  is  just  south  of  Korea, 
and  the  polar  front  is  just  south  of  Kyushu  and  Shikoku.  De¬ 
spite  the  very-  active  weather  systems  the  mean  surface  wind 
speed  over  the  central  Yellow  Sea  in  summer  is  between  5  and 
4  m  s.  which  is  weaker  than  in  winter  (Figure  2b).  June  also 
marks,  historically,  a  jump  in  precipitation  associated  with 
warm,  moist  air  south  of  the  polar  front  [Warn,  1969],  Occa¬ 
sionally.  the  Okhotsk  high  blocks  the  northerly  progression  of 
the  polar  front.  By  July,  however,  high  pressure  (the  Bonin 
high)  to  the  south  and  low  pressure  over  Manchuria  produce 
southerly  winds  carrying  warm,  moist  air  over  the  Yellow  Sea. 
The  summer  monthly  mean  SAT  is  quite  uniform,  around 
24'-26t'C.  and  is  usually  l.5°-2X  warmer  than  the  mean  SST 
\Van  Loon,  1984],  The  wanner  air  causes  a  downward  heat  flux 
at  the  air-ocean  interface.  This  heat  flux  plus  the  strong  down¬ 
ward  net  radiation  stabilizes  the  upper  layer  of  the  water  and 
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Figure  2.  Mean  atmospheric  surface  circulations  in  the  vicinity  of  the  Yellow  Sea  for  (a)  January  and  (b) 
July  (computed  from  the  European  Centre  for  Medium-Range  Weather  Forecasts  (ECMWF)  data). 
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Figure  3.  Eastern  Yellow  Sc:i  (around  36  N)  temperature 
protilcs  during  1950-1988:  (a)  January  and  (b)  June.  Solid 
squares  show  the  observation  stations. 

causes  the  surface  mixed  layer  to  shoal,  creating  a  multilayer 
structure.  Below  ihe  therm ocline.  there  is  a  cold  water  mass, 
commonly  referred  to  as  Yellow  Sea  Bottom  Cold  Water  (YS- 
BCW).  that  remains  unchanged  and  nearly  motionless 
throughout  the  summer  [IJ  and  Yuan.  1992J.  October  is  the 
beginning  of  the  transition  back  to  winter  conditions.  The 
southerly  winds  weaken,  letting  the  sea  surface  slope  reestablish 
the  winter  pattern.  The  SST  steadily  decreases  from  October 
to  January. 


4.  Seasonal  Variation  of  Temperature 
Profiles 

The  Yellow  Sea  shelf  water  masses  experience  a  strong  sea¬ 
sonal  variation  (Figure  3).  Taking  the  eastern  part  of  the  \  el- 
low  Sea  around  36  N  as  an  example,  the  January  historical 
( 1950-1988)  temperature  profiles  (Figure  3a)  show  a  single- 
layer  structure  (i.e..  vertically  uniform  temperature  from  sur¬ 
face  to  bottom).  The  different  lengths  of  these  profiles  in  the 
vertical  are  caused  by  the  different  water  depths  where  the 
observations  were  taken  (see  Figure  I ).  Tins  single-layer  struc¬ 
ture  means  a  very  deep  mixed  layer  extending  from  the  surface 
to  bottom.  On  the  other  hand,  the  June  historical  ( 1950-1988) 
temperature  profiles  (Figure  3b)  indicate  multilayer  structure 
(i.e..  mixed  layer,  thermocline.  and  deep  layer).  There  also 
exists  a  shallow  surface  mixed  layer  in  the  summer  profiles. 
Such  a  strong  seasonal  variation  in  vertical  thermal  structure  is 
caused  by  a  strong  surface  coiling  in  winter  and  a  strong 
surface  warming  in  summer. 

5.  Thermal  Parametric  Model 

During  the  summer  monsoon  season,  most  profiles  in  the 
Yellow  Sea  exhibit  a  mixed  layer,  a  thermocline.  and  a  deep 
layer  (Figure  3b).  which  can  lie  outlined  by  a  "typical  profile. 
To  make  the  model  more  general,  wc  assume  two  deep  layers 
below  the  thermocline  (Figure  4).  When  the  two  deep  layers 
have  the  same  vertical  gradients,  they  become  one  deep  layer. 
If  two  transition  layers  are  added,  the  entrainment  zone  be¬ 
tween  the  mixed  layer  and  thermocline  and  the  transition  zone 
between  the  thermocline  and  the  deep  layer,  the  Yellow  Sea 
thermal  structure  during  summer  can  be  well  resolved.  Wc  use 
a  parametric  model  with  six  layers  [Chu,  1995]  to  diagnose 
shallow-water  multilayer  structure  from  observed  Yellow  Sea 
temperature  profiles.  Each  observed  profile  is  modeled  hv  a  set 
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Figure  4.  The  temperature  and  gradient  space  representations  of  the  features  or  profile  characteristics 
modeled  by  the  feature  model. 
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of  parameters,  most  of  which  have  physical  meaning,  including 
SST.  MLD.  depth  of  the  base  of  the  thermocline.  gradient  in 
the  thermocline  and  deep  layers,  and  additional  parameters 
describing  curvature  between  the  mixed  layer  and  thermocline 
and  curvature  below  the  thermocline.  Among  them.  SST  is 
taken  as  the  observed  values.  The  model  parameters  for  each 
observed  profile  tire  computed  in  gradient  space:  more  specif¬ 
ically.  the  depths  and  gradients  of  the  modeled  features  are  lit 
to  the  vertical  gradient  of  the  observed  profile.  The  parametric 
model  depicts  the  multilayer  structure.  Determination  of  layer 
number  is  based  on  overall  features  of  the  profiles  (Figure  3b). 
In  this  study  the  thermal  parametric  model  consists  of  seven 
depths  (i.e..  six  layers)  and  six  gradients  as  shown  in  Figure  4. 
The  first  and  last  depths  are  assumed  to  lie  at  the  surface  and 
bottom,  respectively,  and  the  gradients  within  the  mixed  layer 
and  thermocline  are  constrained  to  he  constant.  The  gradients 
for  the  four  other  layers  are  assumed  to  vary  with  depth  lin¬ 
early.  The  mean  gradient  is  taken  as  the  representative  value 
for  these  layers.  The  model  parameters  are  calculated  in  the 
gradient  space,  which  will  bring  larger  numerical  errors  due  to 
the  dilfcrentiation.  Wc  will  use  the  optimization  to  filter  out 
the  noise. 


5.1.  Thermal  Parametric  Model 

If  we  consider  profiles  in  the  gradient  space,  i.e.,  G,  •= 
<iT(z)  Hz.  each  profile  can  he  represented  by  the  surface  value 
(SST)  plus  the  gradients,  c.g.. 

[SSI ,  Gf l).  C|L  G;\  C],  r ;).  •  •  •  .  G ft e„_i,  r„)] 


for  the  temperature  profiles.  Here,  ;i  +  1.  is  the  number  of 
data  points,  and  z,(i  =  1.  2.  ••• .  n)  arc  the  depths  of  the 
subsurface  data  points.  For  example.  1(H)  temperature  depth 
points  would  produce  99  gradient  values.  If  the  surface  value  is 
included,  we  have  the  same  amount  of  data  in  the  gradient 
space  as  in  the  original  data  set. 

On  the  basis  of  the  continuity  of  T  and  HT  Hz  at  interfaces 
of  any  two  layers,  a  parametric  model  can  be  constructed  as 

f ""( z)  =  G'i"'z  +  SST  2  S  [  -d,.  0]  ( la) 


(2  +  d,) 
2  idi  -  d,  I 


•l (Gf  +  Gf  %/, -</,)- (Of  -  Gf ) 

•U  +  </:)]+  T’" ■<-</,)  te[-4  -d,l  (1b) 

rS:)  =  Gt' ( 2  +  </;)  +  r™'( -d:>  z  e  [ -d:]  < ic) 

*-•"<*>  =  l'di 1<('T  -  &TU  +  dxst  -  dyGT  i 


temperatures  in  the  mixed  layer,  the  entrainment  zone,  the 
thermocline.  the  transition  zone,  and  the  first  and  the  second 
deep  layers,  respectively.  //  is  the  water  depth,  d,  is  the  mixed- 
layer  depth,  d.  is  the  depth  of  the  thermocline  top.  d-  is  the 
depth  of  the  thermocline  bottom.  d4  is  the  depth  of  the  top  of 
the  first  deep  layer,  and  ds  is  the  bottom  of  the  first  deep  layer 
(Figure  4).  Here  we  assume  constant  vertical  temperature  gra¬ 
dients  in  the  ocean  mixed  layer  (very  small.  Of)  -  0)  and  in 
the  thermocline  (very  large  Gf)  and  linearly  varying  with  - 
in  the  entrainment  zone,  the  transition  zone,  and  the  two  deep 
layers  add  with  average  values  G'f"  Gf.  Gf  \  and  G 
Here  the  mean  gradient  in  the  entrainment  zone  is  the  average 
of  the  mixed-layer  and  thermocline  gradients,  i.e..  Gf'  = 
(G'f”'  T  Gf  )/2.  By  forcing  this  parametric  model  (JaH  If) 
to  each  observed  profile,  we  should  have  a  first  guess  of  the  five 
depths  (d ,.  d..  d ,.  </.,.  d5)  and  a  high  resolution  of  temper¬ 
ature  depth  points  in  the  vertical  in  order  to  obtain  the  five 
temperature  gradients  (G'f”'.  Gf .  Gf,  G'f1'  6^:’). 

Such  a  treatment  provides  the  most  important  features  from 
the  observational  data. 


5.2.  Fligh-Resolution  Profiles  (HP)  Interpolated 
From  Observations 

Each  MOODS  profile  is  linearly  interpolated  to  Ac  =  0.5 
m.  T,  =  Hz.),  where  z,  =2,-,  -  0.5  m  (2,,  =  (J).  The 
N A VOCEANO’s  Digital  Bathymetric  Data  Base  5  (  DBDB5) 
data  set  is  used  to  obtain  water  depth  //.  If  the  five  depths  (d,, 
d2.  d3,  d.,,  d<)  are  known,  we  can  divide  the  data  set  HP  = 
( zr  Tj)  into  six  parts  (mixed  layer,  entrainment  zone,  thermo- 
clinc.  transition  zone,  first  deep  layer,  and  second  deep  layer). 
For  each  layer  we  fit  (  T z;)  to  the  parametric  model  ( I  a)-(  1 0 
and  obtain  a  set  of  temperature  gradients  (Gf.  Gf'',  Gf, 
Gf".  Gf"). 


5.3.  Modeled  Profile  (MP)  Obtained  by  the  Iteration 
Method 

A  modeled  profile  (MP)  with  0.5-m  resolution  can  Ik-  estab¬ 
lished  by  using  the  parametric  model  (la)-(lf)  if  the  five 
depths  (</,.  </,.  d3.  d4.  d. )  arc  given.  In  reality  these  depths 
arc  not  known  prior  to  processing  the  data  and  vary  from  one- 
profile  to  the  other.  Wc  use  the  iteration  method  to  obtain  the 
optimal  MP. 

First,  we  start  with  a  set  of  first-guess  values  of  the  depths 
and  the  five  gradients  (two  constants  and  three  mean  values). 


D""  =  (dr.  d:[,4%,  df,  d'i"), 

Gf  =  <G(r.  of  .  GT.  Gf,  G'f"). 

in  this  studv  we  choose 


(2) 


D  =  (20  m.  24  m.  32  m.  38  m.  H\ 


+  f"*>(-d,)  z  6  [  —dt,  — d3]  (Id) 


G'T  =  tl).  !l.5X/m.  0.005'Cm.  II.  0). 


fJ"(  z)  =  t”'(  -dd  +  (r  +  d,)Gf  -r 
*  [Gr’f-di)  -  Gf  ]  r  e 


(2  -r  d 4  M  2  +•  d,) 
d  5  —  d  i 

l  ~ d„  — dij  (le) 


f"iz)  =  rf-dj  +  (2  +  djGf'  +  -z  +  *s>lj+m 

II  U( 

-  [Gy  1  ( -d<)  -  Gf  ]  :£[-//,  -d5]  (If) 

where  f’”\  P"'\  ff  t'dx\ and  are  modeled 


For  each  HP  profile  we  fit  the  data  (;,.  T, )  to  the  parametric 
model  ( la)— (If)  since  we  know  the  depths  and  the  gradients 
and  obtain  a  zeroth-order  modeled  profile,  called  MP<tM 
( Zj.  fy’1).  The  RMS  error  for  mismatch  of  HP  and  MP1"1  is 
computed  by 


1 

RMS  *  -  «  Y  (,ff  -  Tf 
A  n  ’ 

\  n 

where  k  ().  We  expect  RMS""  to  he  large. 


(3) 
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Figure  5.  Bimonthly  variation  of  SST.  Data  were  obtained  from  the  MOODS  tor  the  years  1950— 1988. 


November 


120  122  124 

Loo^hiiv  tE» 


Second,  we  use  the  iteration  method  to  obtain  optimal  MP 
for  each  HP  profile.  Each  depth  can  only  be  adjusted  one 
vertical  grid  lAr  or  Ac)  for  one  iteration.  From  the  A ih- 
ordcr  l/c  starling  from  l),  the  first  guess)  set  of  depths.  />**’.  we 
have  242(  =  35  —  1)  different  combinations  of  the  depth  ad¬ 
justment. 

D\:  "  =  Da'  +-  RD*\  (4) 

where 

SD',*’=  lAr.  0.  ().  f>.  0). 

8D‘?  =  I  -  Ac.  0.  0.  0.  IT). 

8D%  =  (0.  0.  0.  0.  -Ac). 

We  use  (la)  (If)  to  obtain  242  modeled  profiles,  among 
which  we  pick  a  profile  with  minimum  RMS  error  as  the 
(k  +  I  )th-order  set  of  depths.  Z)**  * 1  \  Wc  repeat  this  pro¬ 
cedure  until  the  minimum  RMS  error  is  achieved.  We  have  two 
check  points  to  terminate  the  iteration:  the  maximum  number 
of  iterations  kmM  and  the  RMS  error  criterion  A, .  At  each 
iteration  (k  <  knMS).  RMS***  is  compared  to  a  user-specified 
criterion  A'. .  If  RMS"1' 1  <  A\  .  we  terminate  the  iteration  and 
obtain  an  optimal  set  of  depths.  If  RMS'*  ’  >  A. .  we  continue 


ihe  iteration  until  k  =  Arnuv.  If  the  RMS  error  at  the  kmx 
iteration  is  still  greater  than  R, .  we  should  reject  the  paramet¬ 
ric  model  ( la)-(lf):  that  is.  the  HP  profile  cannot  lie  fitted  by 
the  parametric  model.  The  rejected  profiles  arc  discarded.  In 
this  study  we  chose  400.  A,  -  0.4°C.  This  R.MS 

criterion  (0.4'C)  was  chosen  due  to  the  accuracy  of  the  tem¬ 
perature  depth  ( c  ().2°C  c2  mi  measured  bv  the  bathythermo¬ 
graph  (BT).  If  the  data  are  obtained  by  more  accurate  instru¬ 
ments  (e.g..  thermometer),  this  criterion  can  be  greatly 
reduced.  The  number  of  rejected  profiles  is  1232.  The  six 
temperature  gradients  were  updated  at  each  iteration.  A  set  of 
the  six  optimal  gradients  is  obtained  when  the  RMS  error  of 
the  temperature  profile  is  less  than  the  criterion. 

6.  Yellow  Sea  Shelf  Thermal  Features 

The  parametric  model  ( I  a )-(  1  f)  transforms  any  profile  (if 
not  rejected)  into  a  set  of  five  depths  (d..  d7.  </  ,.  </  ..  </s  i  and 
six  gradients  O'/"1-  <;</”.  6f*K 

These  parameters  represent  the  most  important  physical  Ica- 
tures  and  are  insensitive  to  the  first-guess  values.  The  first- 
guess  values  only  affect  the  number  of  iterations.  A  good  se¬ 
lection  iff  first-guess  values  reduces  the  number  of  iterations. 
Since  winter  (Deccrnber-April)  profiles  reveal  a  single-layer 
structure  | Fralick.  IW4),  the  set  of  characteristic  parameters 
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Figure  6.  Monthly  variation  of  the  central  Yellow  Sea  (123°- 
124  E.  30  -37'N)  thermal  features:  (a)  SST  (degrees  Celsius) 
and  (b)  MLD  (meters). 

reduces  to  (SST.  H).  Therefore  the  thermal  parametric  model 
was  used  to  process  only  summer  monsoon  season  (May- 
November)  profiles.  Due  to  the  data  sparsity,  the  summer 
profiles  with  the  same  month  and  different  years  were  grouped 
together  to  obtain  the  characteristic  parameters  of  the  monthly 
variations. 

6.1.  Seasonal  Variation  of  Sea  Surface  Temperature 

Bimonthly  evolution  of  SST  (Figure  5).  constructed  from  the 
MOODS  data  (1950-1988).  shows  mime  interesting  features 
which  are  summarized  as  follows. 

1.  During  the  winter  monsoon  season,  SST  has  a  weak 
temporal  variation.  In  the  central  Yellow  Sea.  SST  decreases 
slowly  cooling  rate  ~-4°C/month)  from  November  through 
January,  stays  almost  the  same  ((T-8°C)  from  January  through 
March,  and  increases  slowly  (warming  rate  -3°C'month) 
from  March  through  May.  The  heat  loss  (or  gain)  at  the  ocean 
surface  cools  (or  warms)  the  whole  column  of  the  water,  caus¬ 
ing  a  slow  change  of  SS  T. 

2.  (Tie  winter  SST  (January.  March)  shows  a  tongue-like 
distribution  which  extends  from  southeast  to  northwest,  w  hich 
coincides  with  observations  collected  during  a  single  year 
( 1986)  in  the  Yellow  Sea  \Cireii  etal..  1994],  Such  a  tongue-like 
structure  implies  the  existence  of  the  Yellow  Sea  Warm  Cur¬ 
rent  (YSWC)  in  winter  [Hsueh.  1988). 

3.  During  the  summer  monsoon  season,  on  the  other  hand. 
SS  T  has  a  strong  temporal  variation.  In  the  central  Yellow  Sea. 
SST  increases  rapidly  from  May  through  July  (warming  rate 

?  Cmonth)  and  stays  warm  (above  20  C  )  through  Septem¬ 
ber.  SST  noticeably  cools  from  September  through  December 
(  l  C  month).  The  rapid  change  is  due  to  a  shallower  mixed 


layer  during  the  summer  monsoon  season.  The  heat  gain  (or 
loss)  at  the  ocean  surface  warms  (or  cools)  onlv  the  mixed 
layer,  causing  a  rapid  change  of  SST. 

4.  The  seasonal  variation  of  sea  surface  temperature  has  its 
the  largest  value  (  20X7)  in  the  northern  Yellow  Sea  (north  of 
37eN)  and  decreases  southward. 

Lhe  SSI  over  the  Yellow  Sea  has  a  strong  annual  cycle.  In 
the  central  Yellow  Sea  (I23°-I24°E.  36  -37°N)  the  maximum 
SST  (25.23C)  appears  in  late  July  and  early  August,  and  the 
minimum  SST  (7°C)  occurs  in  March  (Figure  6a). 

6.2.  Mixed-Diver  Depth 

By  late  winter  the  water  throughout  the  Yellow  Sea  is  very 
isothermal.  During  onset  of  the  summer  monsoon  an  increase 
in  the  solar  insolation  and  weakening  winds  begin  the  stratifi¬ 
cation  process  on  the  Yellow  Sea  Shelf.  This  is  supported  by  a 
strong  annual  cycle  of  MLD  in  tite  central  Yellow  Sea  ( 1 23'- 
124  'E.  36  -37°N  ).  During  the  winter  monsoon  season  (Decem¬ 
ber- March).  MLD  has  the  same  value  < — 50  m)  as  the  water 
depth  (single-layer  structure).  During  the  onset  of  the  summer 
monsoon.  MLD  reduces  drastically  from  5(1  to  8  m  (multilayer 
structure  due  to  stratification)  and  then  keeps  small  values 
(less  than  S  m)  throughout  the  whole  summer  season  (Figure 
fib). 

The  Yellow  Sea  bimonthly  progression  of  the  MLD  plots 
(Figure  7)  reveals  some  interesting  features: 

1.  There  are  two  stages  for  the  time  evolution  of  mixed- 
layer  depth:  a  quasi-steady  (and  shallow)  stage  from  May 
through  July  and  an  evident  deepening  stage  from  July  through 
November. 

2.  1  he  contours  of  the  mixed-layer  depth  almost  follow  the 
bathymetry  (Figure  l).The  central  Yellow  Sea  has  a  relatively 
deep  mixed  layer  during  the  whole  summer  monsoon  season, 
whose  depth  is  from  5  to  1(1  m  for  the  quasi-steady  stage  ( under 
very  weak  surface  forcing)  and  deepens  to  20-25  m  during  an 
evident  deepening  stage  caused  by  reduced  insolation  and  in¬ 
creased  wind. 

3.  The  mixed-layer  depth  is  very  shallow  (-5  rn)  in  the 
north  Yellow  Sea  (north  of  37°N)  during  the  quasi-steady  stage 
and  increases  to  15-20  m  during  the  evident  deepening  stage. 
The  shallow  mixed  layer  in  the  first  stage  is  very  sensitive  to  the 
net  surface  heat  flux. 

4.  Dirge  mixed-layer  depth  occurs  during  the  transition 
month  of  November,  when  solar  insolation  has  decreased 
markedly  and  winds  have  picked  up. 

5.  The  mixed-layer  depth  is  also  affected  by  the  large  dis¬ 
charge  from  surrounding  rivers  (Yangtze  River.  Muangho 
River,  etc.).  As  reported  by  Yang  ft  al.  [1983],  the  discharge 
from  the  Yangtze  River  is  maximum  in  July  (16%  of  annual 
discharge)  and  minimum  in  January  (3T  of  annual  discharge). 
Near  the  Chinese  coast  the  mixed-layer  depth  is  minimum  in 
July  (less  than  2  m). 

6.  1  he  mixed-layer  depth  is  strongly  influenced  hv  tidal 
mixing  in  shallow  regions  such  as  along  the  Chinese  coast  in 
water  less  than  20  m  deep.  The  amplitudes  of  M,  and  S,  tides 
there  are  larger  south  of  34°N  than  north  of  >4  \  [Yanagi  and 
Inoue.  1994],  which  implies  stronger  tidal  vertical  mixing  (in 
turn  thicker  mixed  layer)  south  of  34  N.  Hi  is  can  he  clearly 
seen  in  September  and  November  MLD  plots  (Figure  7). 

6-3.  Prominence  of  the  Seasonal  Thermoeline 

July  is  a  month  ol  strong  surface  heating  as  witnessed  by  the 
sharp  increase  in  SST  shown  in  Figure  5.  As  a  result  the 
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Figure  7.  Bimomhly  evolution  of  MLD  (meters)  in  the  Yellow  Sea.  Data  were  gridded  and  horizontally 
interpolated  from  the  MOODS  for  the  years  1950-1988. 


prominence  of  the  thcrmocline  is  expected.  With  the  presence 
of  YSBCW  in  the  central  region  a  positive  temperature  gradi¬ 
ent  across  the  thcrmocline  should  be  apparent  and  appears  in 
the  central  part  of  the  Yellow  Sea  shelf.  The  presence  of  this 
centrally  manifested  gradient  is  more  obvious  in  August  and 
shows  a  stronger  gradient  as  summer  progresses.  September 
shows  a  ridge  in  the  temperature  gradient  distribution  oriented 
north-south  in  the  central  region  which  suggests  the  deeper, 
colder  water  serves  to  maintain  the  temperature  difference 
even  into  the  later  summer.  As  solar  insolation  decreases  in 
October  and  wind  speeds  begin  to  increase,  entrainment  be¬ 
gins  at  the  mixed-layer  base  to  break  down  the  thermocline. 
This  is  borne  out  by  the  October  map  where  the  north-south 
oriented  ridge  in  temperature  gradient  is  still  present  but  the 
magnitude  of  the  gradient  has  decreased  significantly. 

The  temperature  difference  across  the  thermocline.  \Ti'h> , 
is  computed  from  the  output  of  the  parametric  model.  d2.  d^, 
and  G("’\ 

A7,rt'=  G^idi-dJ  (5) 

The  bimonthly  evolution  of  A T<,/"  field  (Figure  8)  indicates 
some  important  features. 

1.  The  contour  of  AT"'*  almost  follows  the  bottom  topog¬ 
raphy  (Figure  1).  No  matter  which  month  in  the  year,  the 


maximum  AT’’''0  always  appears  in  the  central  Yellow  Sea. 
This  phenomenon  can  be  explained  by  the  tidal  vertical  mixing 
[Yanagi  and  Inoue,  1994;  Takahashi  and  Yanae.  1995].  The 
circulations  in  the  Yellow  Sea  during  summer  are  mainly  in¬ 
duced  by  the  sea  surface  heating  and  arc  affected  by  the  ver¬ 
tical  mixing  of  the  tidal  current  and  the  bottom  topography. 
The  homogeneous  water  which  is  formed  by  the  sea  surface 
cooling  in  winter  begins  to  stratify  due  to  the  sea  surface 
heating  in  spring.  However,  due  to  the  horizontal  difference  of 
the  tidal  vertical  mixing  effect  (stronger  tidal  mixing  in  shallow 
shelves),  the  intensity  of  the  stratification  in  the  central  Yellow 
Sea  becomes  stronger  than  the  surrounding  part.  Besides,  sur¬ 
face  heating  during  the  summer  monsoon  season  also  shoals 
the  mixed  layer  and  makes  the  mixed-layer  temperature  (i.e., 
SST)  relatively  uniform  horizontally.  Since  the  minimum  tem¬ 
perature  in  the  deeper  layer  appears  in  the  central  Yellow  Sea 
(i.e.,  YSBCW).  a  strong  thermocline  should  be  located  there. 

2.  The  thermocline  strengthens  during  the  first  half  of  the 
summer  monsoon  season  (May-July).  The  maximum  A 
increases  from  5°C  in  May  to  14°C  in  July.  During  the  second 
half  of  the  summer  monsoon  season  (July-Septcmber)  the 
thermocline  weakens  except  in  the  central  Yellow  Sea  where 
the  thermocline  strength  remains  high  value  (maximum 
A V'h'1  —  14°C).  The  strength  of  the  thermocline  further 


10306 


CHU  ET  .-VL.:  A  PARAMETRIC  MODEL  FOR  YELLOW  SEA  THERMAL  VARIABILITY 


May 


118  120  122  124  126 

Longitude  (E) 


September 


July 


November 


118  120  122  124  126 

Longitude  (E) 


Figure  8.  Thermoclinc  temperature  difference  (degrees  Celsius)  in  the  Yellow  Sea  from  Mav  to  November. 
Data  were  obtained  from  the  MOODS  for  the  years  1950-1988. 


reduces  from  September  through  November  (maximum 
\Toh)  —  5°C).  During  the  winter  monsoon  season  (Decem- 
ber-April)  the  Yellow  Sea  waters  are  well  mixed  vertically,  and 
there  is  no  thermocline  during  that  period. 

7.  Conclusion 

1.  The  thermal  parametric  model  depicted  in  this  paper 
demonstrates  a  good  capability  to  compress  a  large  data  set 
into  a  small  enough  volume  of  data  such  that  the  investigators 
may  readily  assimilate  and  interpret  it.  The  model  clearly 
shows  large  thermal  variabilities  and  the  dual  structure  of  the 
Yellow  Sea  shelf:  a  one-layer  structure  during  the  winter  mon¬ 
soon  season  (December-April)  and  a  multilayer  structure  dur¬ 
ing  the  summer  monsoon  season  (May-Novcmber). 

2.  The  Yellow  Sea  surface  mixed  layer  has  more  thermal 
inertia  during  winter  monsoon  season  than  summer  monsoon 
season.  SST  stays  almost  the  same  from  January-  through 
March  (winter  monsoon)  and,  increases  dramatically  from 
May  through  July  (summer  monsoon  season).  The  maximum 
warming  rate  (7°C/month)  in  summer  is  nearly  2  times  higher 
than  the  cooling  rate  (<  -4°C/month)  in  winter.  'Die  seasonal 
variation  of  SST  has  its  largest  value  (20°C)  in  the  northern 
Yellow  Sea  and  decreases  southward. 

3.  During  the  winter  monsoon  season  the  surface  mixed 
layer  extends  from  the  surface  to  the  bottom.  After  the  onset 


of  the  summer  monsoon  in  May  the  ocean  mixed  layer  shoals 
and  experiences  two  stages  of  evolution.  A  quasi-steady  and 
shallow  stage  (MLD  —  10  m  at  the  central  Yellow  Sea)  appears 
from  May  through  August,  and  an  evident  deepening  stage 
(MLD  ~  20-25  m  at  the  central  Yellow  Sea)  occurs  from 
August  through  November.  During  the  first  stage  the  mixed 
layer  has  the  least  thermal  inertia  and  is  sensitive  to  surface 
thermal  forcing. 

4.  A  strong  thermocline  occurs  during  the  summer  mon¬ 

soon  season,  with  the  largest  \T',h)  occurring  in  September. 
The  thermocline  strengthens  continuously  after  the  onset  of 
the  summer  monsoon  season.  The  maximum  AT"1’1  for  the 
Yellow  Sea  increases  from  5°C  in  May  to  14CC  in  July  and  is 
maintained  (maximum  AT'"0  -  14°C)  through  September. 
The  strength  of  the  thermocline  decreases  from  September 
through  November  (maximum  —  5°C).  During  the 

winter  monsoon  season  the  Yellow  Sea  waters  are  vertically 
well  mixed,  with  little  or  no  stratification. 

5.  The  thermal  structure  is  affected  by  the  air-sea  surface 
heat  flux,  the  wind,  the  tidal  vertical  mixing,  and  the  large 
discharge  from  the  mainland  of  China  (Yangtze  River. 
Huangho  River,  etc.).  The  Yellow  Sea  lies  within  the  Asian 
monsoon  circulation  and  experiences  hot.  humid  summers  and 
long,  cold,  dry  winters.  In  winter  the  Yellow  Sea  is  dominated 
by  the  cold  northerly  wind  associated  with  the  Siberian  high. 
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The  SST  is  6'C  ai  the  northern  extent  and  1(YC  at  the  south¬ 
east  extent.  The  January  to  surface  air  temperature  (SAT) 
varies  from  0°  to  8°C  in  the  Yellow  Sea  [ Chu  et  al..  1997]. 
roughly  2°-6°C  cooler  than  SST.  The  Yellow  Sea  surface  loses 
tremendous  heat  to  the  atmosphere.  The  upward  buoyancy 
flux  at  the  air-ocean  interface  (  thermal  forcing),  together  with 
the  strong  wind  stress  (mechanical  forcing),  generates  turbu¬ 
lence  and  mixes  the  surface  water  with  the  deeper  water.  The 
mixed  layer  is  at  its  deepest  (it  usually  fills  the  whole  water 
column)  during  winter  owing  to  both  convection  and  wind 
mixing  by  the  strong  northeast  monsoon  winds.  In  summer  the 
Yellow  Sea  is  dominated  by  the  warm  and  moist  southwesterly 
wind  associated  with  the  atmospheric  low-pressure  system  over 
Asia  (the  Manchurian  low)  and  high-pressure  system  in  the 
southeast  (the  Bonin  high).  The  summer  SAT  is  quite  uniform 
and  around  24*-26'C  [see  Chu  et  al ..  1997].  It  is  usually  1.5*— 
2*C  warmer  than  SST  [Van  Loon.  1984].  The  warmer  air  brings 
downward  heat  flux  at  the  air-ocean  interface.  This  heal  flux 
plus  the  strong  downward  net  radiation  stabilizes  the  upper 
layer  of  the  water  and  causes  the  surface  mixed  layer  to  shoal. 
Furthermore,  near  the  Chinese  coast  the  MLD  has  a  minimum 
in  July  (less  than  2  m).  This  is  mainly  due  to  the  maximum 
freshwater  influx  at  the  major  river  mouths  in  July.  The  mixed- 
layer  depth  at  the  southwest  comer  of  the  Yellow  Sea  (near  the 
Yangtze  River  mouth)  increases  from  2  m  in  July  to  the  whole 
water  column  in  November. 

6.  This  parametric  model  can  also  be  used  to  process  out¬ 
put  from  any  dynamical  model  and  to  obtain  meaningful  prod¬ 
ucts.  i.e..  MLD.  thermocline  depth,  thermocline  strength,  and 
deep  layer  stratification. 
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